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1 Notations and preliminary remarks 

Selfsimilarity is a very popular research topics since a few years, and it has been ap- 
proached from many different standpoints producing an unavoidable level of confu- 
sion [1]. The present paper is devoted to a short summary of the properties of some 
important and well known families of laws: infinitely divisible, selfdecomposable and 
stable (for details see for example [H El El [5] ) • First of all we will recall their role in 
the formulation and in the solutions of the central limit problem: as we will see in 
the next section this amounts to a quest for all the limit laws of sums of independent 
random variables. Our families of distributions will then be analyzed by means of both 
their possible decompositions in other laws, and the explicit form of their characteristic 
functions: the celebrated Levy-Khintchin formula. In particular it will be discussed 
the intermediate role played by the selfdecomposable laws between the more popular 
stable, and infinitely divisible distributions. We will then explore these laws in con- 
nection with the additive and the Levy processes, looking for their importance with 
respect to both the properties of stationarity and selfsimilarity. In particular it will be 
recalled how from selfdecomposable distributions it is always possible to define both 
stationary and selfsimilar additive processes which - with the exception of important 
particular cases - will in general be different. A few remarks are also added to show the 
differences between this selfsimilarity and that of the well known fractional Brownian 
motion. We will also analyze the Ornstein-Uhlenbeck processes driven by these Levy 
noises, and their stationary distributions which are always selfdecomposable. We will 
finally elaborate a few examples to illustrate these results, and to compare the behavior 
of our processes. 

In what follows we will adopt the following notations (for details see for exam- 
ple [21 16]): X, Y, . . . will denote the random variables, namely the measurable functions 
X(lo),Y(u), . . . defined on a probability space (S7, J 7 , P) where f2 is a sample space, T 
is a iT-algebra of events and P a probability measure. We will then respectively write 
F(x), G(y), ... for their cumulative distribution functions 



where the symbol E denotes the expectation value of a random variable according to 
the probability P, namely for example 



When they exist, f(x),g(y), . . . will be the probability density functions, and in that 
event we will have 



F{x) = P{X < x} , G(y) = P{Y < y} 



and (f(u), x( v )i ■ ■ ■ f° r their characteristic functions 




x F(dx) . 






A law will be indifferently specified either by its cumulative distribution function (or 
density function), or by its characteristic function. To say that a random variable X is 
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distributed according to a given law we will also write either X ~ F(x) or X ~ f(u)', 
if a family of laws is denoted by a specific symbol, say C, then we also write X ~ C to 
say that X is distributed according to one of these laws. When two random variables 
X and Y are identically distributed we will adopt the notation 

X = Y . 

If a property is true with probability 1 we will also say that it is true P-almost surely 
and we will adopt the notation P-a.s. A type of laws (see [2] Section 14) is a family 
of laws that only differ by a centering and a rescaling: in other words, if f(u) is 
the characteristic function of a law, all the laws of the same type have characteristic 
functions e lbu (f(au) with a centering parameter b G R, and a scaling parameter a > 
(we exclude here the sign inversions). In terms of random variables this means that 
the laws of X and aX + b (for a > 0, and t £ 1) always are of the same type, and 
on the other hand that X and Y belong to the same type if it is possible to find 
a > 0, and i G K such that Y = aX + b. We will also say that a random variable 
X ~ <p(u) and its law are composed of X\ ~ <fi(u) and X<i ~ ^(u) when X\ and X2 

are independent and X = X\ + X2, or equivalently when ip(u) = (pi(u)(p2(u); then 
X\ and X2 are also called components of X. Of course, in terms of distributions, a 
composition amounts to a convolution; then, for example, if Af(a 2 ,b) denotes a normal 
law with expectation b and variance a 2 , the composition of two normal laws will also be 
indicated as M{a\,b\) * A^a^,^)- For the stochastic processes X(t),Y(t), ... we will 
say that X(t) and Y(t) are identical in law when the systems of their finite-dimensional 
distributions are identical, and in this case we will write 

X(t) = Y{t) . 

A process X(t) is said to be stochastically continuous if for every e > 
jim o P{|X(t + At) -X(t)\ > e} = 0, Vt > 0. 

where P{. . .} denotes the probability for the increment \X(t + At) — X(t)\ of being 
larger that e > 0. In this paper all our random variables and processes will be one 
dimensional. 

In this exposition we do not pretend neither rigor, nor completeness: we just list the 
results and the properties that are important to compare, and we refer to the existing 
literature for proofs and details, and for a few hints about possible recent applications. 
In fact the aim of this paper is just to draw an outline showing - without embarrassing 
the reader with excessive technical details - the deep, but otherwise simple ideas which 
are behind the properties of our families of laws and the simplest procedures to build 
from them the most important classes of processes. We hope that this, with the aid 
of some telling examples, will be helpful to approach this field of research by clarifying 
the roles, the differences and the subtleties of these laws and processes. 
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2 The Central Limit Problem 

2.1 The classical limit theorems 

It is well known that there are three kinds of classical limit theorems (all along this 
paper wherever we speak of convergence it is understood that we speak of convergence 
in law) which are characterized by the form of the respective limit laws: 

• the Law of Large Numbers with limit laws of the degenerate type 5^ with charac- 
teristic function 

(p(u) = e ibu ; 

• the Normal Central Limit Theorem whose limit laws are of the gaussian type 
M(a 2 ,b) with 

<p(u) = e ibu e -a 2 u 2 /2. 

• and the Poisson Theorem whose limit laws are of the Poisson types V(X;a,b) 
with 

^{u) = e lbu e x ^ au ' l l (1) 

The exact statements of these theorems in their traditional formulations are reprinted 
in every handbook of probability (see for example [6] p. 323 and following), and we 
will not reproduce them once more. Remark however that in our list, while speaking 
of the type for the degenerate and the gaussian laws, we also referred to the types for 
the Poisson laws. In fact all the normal laws M(a 2 ,b) with expectation 6 6 R and 
variance a 2 > constitute a unique type, and the same is true for the family of all 
the degenerate laws. On the other hand the standard Poisson laws V(X) = V(X; 1,0) 
with different parameters A belong to different types: indeed we can not recover a law 
V(X) from another V(X') (with A ^ A') just by means of a centering and a rescaling. 
In other words we could say that every Poisson law V(X) with a given A generates - by 
centering and rescaling - a distinct type V(X;a,b) with characteristic functions (pj. 

The particularities of the Poisson laws with respect to the other two families of limit 
laws are also apparent from their properties of composition and decomposition. For a 
family of laws (not necessarily a type) to be closed under composition means that the 
composition (convolution) of two laws of that family still belongs to the same family. 
On the other hand closure under decomposition means that if a law of the family is 
decomposed in two laws, these two components necessarily belong to the same family. 
The types of limit laws appearing in the classical limit theorems show an important 
form of closure under composition and decomposition summarized in the following 
result (see [2] p. 283): the degenerate and normal types are closed under compositions 
and under decompositions; the same is true for every family of Poisson laws V(X;a,b) 
with the same a. The closure under composition is in fact an elementary property; not 
so for the closure under decomposition: the proofs for the normal and the Poisson case 
were given in 1935-37 by H. Cramer and D.A. Raikov respectively. The normal and 
Poisson composition and decomposition properties can then be stated by saying that 

N(a\M) *N{a 2 2 ,b 2 ) = M{a 2 l + a 2 2 ,b l + b2) 
r(X 1 ;a,b 1 )*r(X 2 ;a,b 2 ) = 7>(Ai + A 2 ; a, h + b 2 ) 
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Hence it is also apparent that a type of Poisson laws never is closed under composition 
and decomposition: we are always obliged to switch from a Poisson type to another 
while composing and decomposing them. 

It is important to remark, however, that the classical limit theorems are not the 
embodiment of these composition and decomposition properties only, but they are more 
far-reaching and profound statements. In fact not only these theorems deal with limits 
of sums of independent random variables of a type which in general is different from 
that of the eventual limit laws - the Poisson law is the limit of sums of Bernoulli 0- 
1 random variables, while the normal and degenerate laws are limits of sums of still 
more general random variables - but also the distribution of the sum of the n random 
variables does not coincide with the limit law at every step n of the limiting process, 
as happens instead in a simple decomposition. 

2.2 Formulations of the Central Limit Problem 

To formulate the Central Limit Problem (CLP) let us look at it first of all in terms of 
sequences of random variables: usually we take a sequence of random variables, and 
then the sequence of their sums S n = X\ + . . . + X n . In this case when we go from S n 
to S n +\ we just add another random variable X n without changing the previous sum 
S n . However this is not the more general way to produce sequences of sums of random 

(n) 

variables. Consider indeed a triangular array of random variables X} 

X? 

y(2) y{2) 



y(n) Y (n) y(n) 
yv-L , i ■ ■ ■ ) 

with n = 1,2,... and k = 1, 2, . . . , n, and suppose that 

(n) (n) 

1. in every row n £ N the random variables X\ , . . . , Xn are independent, 

(n) 

2. the X^, are uniformly, asymptotically negligible, namely that 

max P{|^ n) | > e} 0, Ve>0. 

Define now the consecutive sums 

n 

^ = £4^ (2) 
fc=l 

Then the central limit problem for consecutive sums of independent random variables 
(CLP\) reads: find the family of all the limit laws of the consecutive sums ([2]) and the 
corresponding convergence conditions (see [2] p. 301-2). Three remarks are in order 
here: 
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• for a given n S N the random variables X\ , . . . , X„ are independent, but in 
general they are neither identically distributed, nor of the same type; 

• going from the row n to the row n + 1 the random variables and their laws change: 
in general X% and X^ +1 ^ are neither identically distributed, nor of the same 
type; as a consequence going from S n to S n +i we not only add the ra+l-th random 
variable, but we also are obliged to adjourn the law of S n ; 

• the uniformly, asymptotically negligible condition is an important technical re- 
quirement added to avoid trivial answers to the CLP\ (see for example [2] p. 302); 
in fact without this condition it is easy to show that any law cp would be limit law 
of the consecutive sums ([2]): it would be enough for every n to take X^ ~ ip, and 

(n) 

X^, = P-a.s. for k > 1. The uniformly, asymptotically negligible condition 
will tacitly be assumed all along this paper. 

Important particular cases of the CLP\ are then selected when we specialize the se- 
quence Xjf 1 in the following way (see [2] p. 331): let us suppose that there is a sequence 
Xk, k = 1,2,..., of independent (but in general not identically distributed) random 
variables, and two sequences of numbers a n > and b n E M, n = 1, 2, . . ., such that for 
every k and n 

X^ = — (x k -^). (3) 
a n V n J 

It is apparent that now going from n to n' we just get a centering and a rescaling of every 
Xk, so that X^ and X^ 1 ^ always are of the same type. Then the consecutive sums 
take the form of normed sums (namely centered and rescaled sums) of independent 
random variables 

Sn = t 4"' = t f U ~ 'A = ± (t X k - K ) = ^ (4) 

k=l k=l ° n v n/ an \k=\ J an 

where we adopt the notation 

n 

Sn = x^. 

k=l 

Then the central limit problem for normed sums of independent random variables 
(CLP2) reads: find the family of all the limit laws of the normed sums (JH) and the 
corresponding convergence conditions. Finally there is a still more specialized formula- 
tion of the central limit problem when we add the hypothesis that the random variables 
Xfc are not only independent, but also identically distributed (see [2] p. 338): in this 
case we speak of a central limit problem for normed sums of independent and identically 
distributed random variables (CLP3). 

2.3 Solutions of the Central Limit Problem 

The answers to the different formulations of the central limit problem need the definition 
of several important families of laws that are much more general than the Gaussian type, 
and that can be defined by means of the properties of their characteristic functions. 



N Cufaro Petroni: Selfdecomposability and selfsimilarity 



7 



Suppose that (p(u) is the characteristic function of a law: we will say that this law 
is infinitely divisible (see [2j] p. 308) if for every n G N we can always find another 
characteristic function <p n (u) such that 

<p(y) = & n (u)] n . 

Apparently the name comes from the fact that our law can always be decomposed in 
an arbitrary number of identical laws; remark however that for different values of n we 
get in general laws (f n (u) of different types. In terms of random variables if X ~ '■piu) 
is infinitely divisible, then for every n E N we can find n independent and identically 
distributed random variables x[ n \ . . . ,xj% all distributed as <p n {u) and such that 

n 

x±Y,xl"\ 

k=l 

namely X is always decomposable (in distribution) in the sum of an arbitrary, fi- 
nite number of independent and identically distributed random variables. Let us call 
J- id the family of all the infinitely divisible laws. Many important distributions are 
infinitely divisible: degenerate, Gaussian, Poisson, compound Poisson, geometric, Stu- 
dent, Gamma, exponential and Laplace are infinitely divisible. On the other hand the 
uniform, Beta and binomial laws are not infinitely divisible: in fact no distribution 
(other than the degenerate) with bounded support can be infinitely divisible (see [5] 
p. 31). Remark that if <f(u) is an infinitely divisible characteristic function, then also 
ip^(u) is an infinitely divisible characteristic function for every A > (see [5] p. 35): 
we will see that this is instrumental to connect the infinitely divisible laws to the Levy 
processes. 

A second important family of laws selected by their decomposition properties is 
that of the self decomposable laws (see [2] p. 334): a law (p(u) is selfdecomposable when 
for every a £ (0, 1) we can always find another characteristic function (p a {u) such that 

<p(u) = (p(au)<p a (u). (5) 

In terms of random variables this means that if X ~ p(u) is selfdecomposable, then 
for every a £ (0, 1) we can always find two independent random variables, X' = X and 
Y a ~ (faiu), such that 

X = aX' + Y a . 

In other words for every a £ (0, 1) X can always be decomposed into two independent 
random variables such that one of them is of the same type of X. It can be shown that 
every selfdecomposable law, along with all its components, is also infinitely divisible 
(see p. 335), so that if we call !Fsd the family of all the selfdecomposable laws, 
then J~sd ^ J~id- The Gaussian, Student, Gamma, exponential and Laplace laws are 
examples of selfdecomposable laws (see [5] p. 98). On the other hand the Poisson laws 
are not selfdecomposable: they only are infinitely divisible. 

Finally we will say that a law (p(u) is stable (see |2] p. 338) if for every c\ > and 
C2 > we can find a > and b £ R such that 



e l u ip(au) = ip{c\u)ip(c2u). 
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This means now that if X ~ ty?(w) is stable, then for every c\ > and C2 > we can 

find two independent random variables X\ = X and X2 = X, and two number a > 
and 6 such that 

+ 6 = c x Xi + c 2 X 2 . 

In other words we can always decompose a stable law in other laws which are of the 
same type as the initial one. Our definition can also be reformulated in a slightly 
different way (see [5] p. 69): for every c > we can always find a > and b G M such 
that 

e*V(a U ) = [^(u)] c , 

and this apparently means that for every c > the law c is of the same type as 

(p(u). A law is also said strictly stable if for every c > always exists a > such that 

tp(au) = [<p(u)] c . (6) 

All the stable laws are selfdecomposable, and hence if Tst is the family of all the stable 
laws we will have Tst Q 3~SD- in fact our classification of laws in only three families 
(infinitely divisible, selfdecomposable and stable) is an oversimplification of a much 
richer structure explored for example in [5], Chapter 3. Among the classical laws only 
the Gaussian and the Cauchy laws are stable. 

Remark that if a law <p belongs to one of our families, then also all its type belongs 
to the same family. Hence it would be more suitable to say that Tip, ^sd and Tst 
are families of types of laws; however in the following, for the sake of simplicity, this 
will be understood without saying. The relevance of our three families of laws lies in 
the fact that they represent the answers to the three formulations of the central limit 
problem discussed in the previous section. In fact it can be shown that Tib , Tsd and 
Tst exactly coincide with the families of the limit laws sought for respectively in CLP\, 
CLP2 and CLP3 (see for example [2] p. 321, p. 335 and p. 339 for the three statements). 
To summarize these results - a few examples will be shown in the Section 14.21 - we can 
then say that: 

• the laws of the normed sums ([!]) of sequences Xj~ of independent and identically 
distributed random variables converge toward stable laws; in particular: when 
the Xk have finite variance their normed sums converge - according to the clas- 
sical theorem - to normal laws, while the non Gaussian, stable distributions are 
limit laws only for sums of random variables with infinite variance; a well known 
example of the non Gaussian limit laws is the Cauchy distribution; 

• the laws of the normed sums (J4]) of sequences X^ of independent, but not neces- 
sarily identically distributed, random variables converge toward selfdecomposable 
laws; the special case of the stable laws is recovered when the X^ are also identi- 
cally distributed; in other words when a selfdecomposable law is not stable it can 
not be the limit law of normed sums of independent and identically distributed 
random variables; 

• finally the laws of the consecutive sums ([2]) of triangular arrays X^ converge to 
infinitely divisible laws: a classical example of non selfdecomposable convergence 
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is represented by the Poisson Limit Theorem recalled in the Section I2.lt °f course 
the selfdecomposable case is obtained when the triangular array has the form ([3]) 
and the consecutive sums become normed sums (j3J); however infinitely divisible 
laws which are not selfdecomposable (as the Poisson law) can not be limit laws 
of normed sums of independent random variables. 

Nothing forbids, of course, that in this scheme a Gaussian law be also the limit law for 
consecutive sums of triangular arrays that do not reduces to normed sums. To complete 
the picture we will hence just recall here that there are also necessary and sufficient 
conditions for the convergence to normal laws of triangular arrays of independent, but 
not necessarily identically distributed, random variables with finite variances (see [6] 



Since the stable distributions are limit laws of normed sums (|4]) of independent and 
identically distributed random variables we can also introduce the notion of domain 
of attraction of a stable law (pg : we will say that a law ip belongs to the domain of 
attraction of ips when we can find two sequences of numbers, a n > and b n , such that 
the normed sums @ of a sequence of random variables all distributed as (p, converge 
to ipg. Remark that this definition can not be immediately extended to the non stable 
distributions which are not limit laws of normed sums of independent and identically 
distributed random variables, so that we can not speak of a unique distribution ip being 
attracted by the limit law. 

Every law with finite variance belongs to the domain of attraction of the normal 
law (see [2] p. 363). It is also important to stress here that, while all stable laws are 
attracted by themselves (see [2] p. 363), a non stable, infinitely divisible law (p - which 
always is in itself the limit law of a suitable consecutive sum ([2]) of some triangular 
array of random variables - also belongs to the domain of attraction of some stable 
law (pg: normed sums of random variables all distributed as <p will converge toward 
some stable law ipg. For instance the Poisson law - which is an infinitely divisible limit 
law, as the Poisson Theorem shows - apparently also is in the domain of attraction of 
the normal law since it has a finite variance: a normed sum of random variables all 
distributed according to the same Poisson law will converge to the Gauss law. Finally 
we recall, without going into more detail, that for a given law ip it is always possible 
to find if it belongs to some domain of attraction, and then it is also possible to find 
both the stable limit law (pg, and the admissible numerical sequences a n > and b n 
entering in the normed sums (j4|) converging to (pg (see [2] p. 364). 

2.4 The Levy-Khintchin formula 

It is not easy to find out if a given law (p(u) belongs to one of the families defined in the 
previous section just by looking at the definitions introduced up to now. It is important 
then to recall the explicit form of the characteristic functions of our families of laws 
given by the Levy-Khintchin formula. It can be proved (see [2J p. 343) indeed that 
the logarithmic characteristic ip(u) = log ip(u) of an infinitely divisible law is uniquely 
associated, through the formula 



p. 326). 
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to a generating triplet (/3 2 , L, 7) where 7, (3 € R, and the Xewy function L(x) is defined 
on R\{0}, is non decreasing on (— 00, 0) and (0, +00), with L(±oo) = and 



lim / y 2 dL(y) < +00 

Js<\v\<x 



JS<\y\<x 

for x > finite. In other words a law will be infinitely divisible if and only if its 
logarithmic characteristic satisfies the relation ([7]) for a suitable choice of (/? 2 ,L,7). 
The Levy function L(x) also defines the Levy measure (for details see [5] Section 8) 



v(B) = [ dL{x) 
Jb 



for every measurable set B of R and v ({0}) = 0, so that often we will refer to (/3 2 , v, 7) 
as the generating triplet. When the Levy measure v is absolutely continuous we will 
also denote by W(x) = L'(x) its density. Remark that the Levy-Khintchin formula (J7J) 
can be given in a variety of equivalent versions (see [5] p. 37) by suitably choosing the 
integrand functions, but we will not go into such details. 

For a given law even the verification of the formula (|7|) is not in general an easy task. 
In the case of stable laws, however, the Levy-Khintchin formula is considerably simpler 
since it no longer explicitely involves integrals on the Levy measure. A stable law is 
characterized by a parameter < a < 2 (see [5] p. 76) and it is then also said a-stable 
(a-stable): the case a = 2 corresponds to the Gaussian laws, and to a vanishing Levy 
measure. In the non Gaussian a-stable cases (0 < a < 2) on the other hand the Levy 
measure is not zero, it is absolutely continuous and we have (see [5] p. 80) 



dL(x) = W(x) dx 



Ax 1 a dx, for x > 0, 
B\x\~ 1 ~ a dx, for x < 0, 



with A > 0, B > and A + B > 0. Finally in both cases - Gaussian and non Gaussian 
- the logarithmic characteristic must satisfy the following relation (see [5] p. 86) 



iau — b\u\ a (l — i sign (u) ctan if a 7^ 1, 

iau — b\u\ (l + i sign (u) - clog |w|) if a = 1, 



where a G (0, 2], a £ R, b > and |c| < 1. The Gaussian case simply corresponds 
to a = 2. When the law is also symmetric the characteristic function is real and the 
formula reduces itself to the quite elementary expression 

= e ~ b \ u \ a , 0<a<2. (8) 

The form of the Levy-Khintchin formula, or equivalently of the triplet (/3 2 ,L,7), 
of the selfdecomposable laws, on the other hand, is not so simple. They in fact play in 
some sense a sort of intermediate role between the generality of the infinitely divisible 
laws and the special properties of the stable laws. It can be proved indeed (see [5] p. 95) 
that a law is selfdecomposable if and only if its Levy measure is absolutely continuous 
and its density is 

W(x) = L'(x) = ^ (9) 
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where the function k(x) is non negative, is increasing on (— oo,0) and is decreasing 
on (0, +00). By the way this also show why a Poisson law (whose Levy measure is 
not absolutely continuous) can not be selfdecomposable. As a consequence the Levy- 
Khintchin formula (|7|) of the selfdecomposable laws is more specialized than that of 
the general infinitely divisible laws, but it still contains a non elementary integral part. 

3 Levy processes and additive processes 

An additive process X(t) is a stochastically continuous process with independent incre- 
ments and X(0) = , P-a.s. (namely the probability of not being zero vanishes); on the 
other hand a Levy process is an additive process with the further requirement that the 
increments must be stationary (for further details see [5] Section 1). The stationarity of 
the increments means that the law of X(s + 1) — X(s) does not depend on s. Processes 
with independent increments are also Markov processes, and hence the entire family of 
their joint laws at an arbitrary, finite number of times can be deduced just from the 
one- and the two-times distributions. In other words it is enough to know the laws of 
the increments to have the complete law of the process. This of course is a very good 
reason to be interested in Markov, and in particular in additive processes, but it must 
be recalled here that there are also non Markovian processes which can still be defined 
by means of very simple tools. An important example that will be briefly mentioned 
later is the fractional Brownian motion which takes advantage of being a Gaussian 
process to make up for its lack of Markovianity. It is also important to recall here that 
- with the exception of Gaussian processes - the additive processes trajectories can 
make jumps. This does not contradict their stochastic continuity because the jumping 
times are random, and hence, for every t, the probability of a jump occurring exactly 
at t is zero. 

3.1 Stationarity and infinitely divisible laws 

In the following the law of the process increment X(t) — X(s) will be given by means of 
its characteristic function (fr s j(u), so that the stationarity of the Levy processes simply 
entails that <j> s t{u) only depends on the difference r = t — s: in this case we will use 
the shorthand, one-time notation (p T {u). As for every Markov process the laws of the 
increments of an additive process must satisfy the Chapman-Kolmogorov equations 
which for the characteristic functions are 

<t>r,t{u) = (j>r,s{u)(j)s,t{u) , < T < S < t ] (10) 

for Levy (stationary) processes these equations take the form 

4>a+r(u) = (j) c {u)(j) T {u) , CJ,r>0. (11) 

There is now a very simple and intuitive procedure to build a Levy process: take 
the characteristic function (p(u) of a law and define 



<f>t(u) = [tp(u)] 



(12) 
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It is immediate to see that <pt(u) satisfies (jlip . so that it can surely be taken as the 
characteristic function of the stationary increments of a Levy process. Here T plays the 
role of a dimensional time constant (a time scale) introduced to have a dimensionless 
exponent. To have a consistent procedure, however, we must be sure that when <p(u) is 
a characteristic function, also (f>t(u) in (|12p is a characteristic function for every t > 0, 
but unfortunately this is simply not true for every characteristic function (p(u). We are 
led hence to ask for what kind of characteristic functions (|12j) is again a characteristic 
function. We know, on the other hand, that if <p is an infinitely divisible characteristic 
function, then also tp x with A > is a characteristic function, and an infinitely divisible 
one too. In fact it is possible to show that Q12|) is a characteristic function if and only 
if ip(u) is infinitely divisible. In other words there is a one-to-one relation between 
the class Fid of the infinitely divisible laws and that of the Levy processes (see [5] 
Section 7). Remark however that in general for a Levy process defined by (j!2|) the 
infinitely divisible law of the increments at a generic time t is neither <p(u), nor of the 
same type of <p(u). Only at t = T the law is necessarily ip(u), while for t ^ T it can be 
rather different and - but for few well known cases - its explicit cumulative distribution 
function (or density function) could be quite difficult to find. 

On the other hand when cp(u) is a stable law it is easy to see from the very defi- 
nition (jSJ) of stability that at every time t the law of the increments (|12p will always 
belong to the same type (this is famously what happens for the Gauss and Cauchy 
laws). In this case we speak of a stable process, and its evolution can be summarized 
just in the time dependence of the law parameters which will produce a trajectory 
inside a unique type. A different situation arises instead when (p(u) only belongs to a 
family of infinitely divisible laws closed under composition and decomposition. As we 
have already remarked these families do not in general constitute a type (as the family 
of the Poisson laws 'P(A)): if however they are closed under composition (as are both 
the Poisson and the Compound Poisson processes) the law p2|) of the increment of the 
Levy process stays in the same family of laws all along an evolution which is described 
by the time dependence of the law parameters; this however does not amount to the 
stability of the process since our family is not a single type. 

Since every Levy process is associated to an infinitely divisible law <p(u) = e^ u \ and 
since every infinitely divisible law is associated to a generating triplet (/3 2 , v, 7) we will 
also speak of the logarithmic characteristic ip( u ) and of the generating triplet (/3 2 , v, 7) 
of a Levy process. In this case however the Levy measure v has also an important 
probabilistic meaning w.r.t. the Levy process (see for example [7] pp. 75-85): for every 
Borel set A of R ; u(A) represents the expected number, per unit time, of (non-zero) 
jumps with size belonging to A. It can also be proved that for every compact set A 
such that ^ A we have u{A) < +00, namely the number of jumps per unit time of 
finite (neither infinite, nor infinitesimal) size is finite. Remark however that this does 
not mean that v is a finite measure on R: in fact the function L(x) associated to v can 
diverge in x = so that the process can have an infinite number of infinitesimal jumps 
in every compact [0, T\. In this case, when i/(R) = +00, we speak of an infinite activity 
process, and the set of the jump times of every trajectory will be countably infinite and 
dense in [0, +00]. For the sake of simplicity we will not introduce here the important 
Levy-Ito decomposition of a Levy process into its continuous (Gaussian) and jumping 
(Poisson) parts: the readers are referred to [7], Section 3.4 for a synthetic treatment. 
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3.2 Selfsimilarity and selfdecomposable laws 

A process X(t) (possibly neither additive, nor stationary) is said to be self similar when 
for every given a > we can find b > such that 

X{at) = bX(t), 

namely when every change a in the time scale can be compensated in distribution by 
a corresponding change b in the space scale. In terms of the increment characteristic 
functions this means that for every a > we must have a b > such that 

4>as >a t{u) = (f>s,tO>u). (13) 

In fact it can be proved more about the form of this space-time compensation: given 
a selfsimilar process we can always find H > such that b = a H (see [5] p. 73). This 
number H is called the exponent of the process or Hurst index, and we will also speak 
of H -selfsimilar processes. For further details about selfsimilar, additive processes see 
also 0. 

Since a Levy process is completely specified by f 1 1 2 [) as characteristic function of 
its increments, then in this case the selfsimilarity means that for every a > it exists 
b > such that 

[ip(u)] at / T = [<p(bu)f T . (14) 



From the definition (|6|) of the strictly stable laws and from (114ft it is easy to understand 
then that the unique selfsimilar Levy processes must be strictly stable. For instance in 
a Wiener process we have (p(u) = e~~ u a / 2 , namely 

Vp(u)\^ = e~^l\ D = ^, 

and hence 

[ v {u)] at ' T = e ~ u2Dat / 2 

so that b = yfa (namely H = 1/2) is the required compensation. This means that, 
insofar as the coefficient D remains the same, we can change the space and time scales 
a and T (namely we can change the units of measure) without changing the Wiener 
process distribution. More precisely than these simple remarks, it can be proved that 
a Levy process X(t) is selfsimilar if and only if it is strictly stable (see [5] p. 71). 

Things are rather different, however, when we consider only additive (not necessarily 
Levy) processes, namely when we can also live without stationarity. Now we must stick 
to the general selfsimilarity equation (|13p . and we must remark again that there is 
another simple, intuitive procedure producing additive (but not necessarily stationary), 
selfsimilar processes: simply consider a characteristic function tp(u), a real number 
H > and take 

M«) = ~7 Tt4- (15) 



It is now apparent that the characteristic functions of this family satisfy the equa- 
tion (jlUp and are also selfsimilar according to the definition (fT3"j) with the space-time 



N Cufaro Petroni: Selfdecomposability and selfsimilarity 



14 



scale compensation produced by b = a . Namely f)15|) produces an .H-selfsimilar pro- 
cess. Of course we must ask here the same question surfaced w.r.t. equation (|12p in 
the case of stationary processes: when can we be sure that the function defined by 
the ratio (|15p of two characteristic function still is the bona fide characteristic function 
of a law? Even in this case, however, the answer can be hinted to by looking at the 
definition ([5|) of a selfdecomposable characteristic function . More precisely it can be 
proved that (|15p is a characteristic function if and only if (p(u) is selfdecomposable 
(see [5] p. 99): this ultimately brings out the intimate relation connecting selfsimilarity 
and selfdecomposability. 

Since selfdecomposable laws are also infinitely divisible the previous remarks show 
that from a given selfdecomposable <p(u) we can always produce two different kinds 
of processes: a Levy process whose stationary increments follow the law <\12\) ; and a 
family of additive, selfsimilar process - one for every value of H > - whose (possibly 
non stationary) increments follow the law (|15p . All these processes generated from the 
same <p(u) are in general different with one exception: when <p(u) is an a-stable law the 
associated Levy process coincide with the if -selfsimilar one with Hurst index H = 1/a. 
In this last case indeed the characteristic functions (|12p and (|15p are identical (to see it 
take for example the symmetric form ([8]) of a stable characteristic function). Remark 
also that the index of an a-stable law always satisfies 0<a<2(a = 2 for the Gaussian 
law) , and that this is coherent with the limitation H = 1/a > 1/2 for the Hurst index 
of the stable, selfsimilar processes (see [5] p. 75). On the other hand, in every other 
case (either non-stable, or a-stable with a ^ \/H), from a selfdecomposable law <p(u) 
we can always build a Levy, non selfsimilar process from (|12p . and a family of additive, 
selfsimilar processes with non stationary increments from (|15p . Finally from a infinitely 
divisible, but not selfdecomposable law we can only get a Levy process from (fT2j) , but no 
selfsimilarity is allowed. For more details about present interest of the selfdecomposable 
distributions and selfsimilar processes in the applications see for instance [9] and |10j . 

Remark that selfsimilarity is not tied to the dependence or independence of the 
increments: we have seen here that among independent increment processes we find 
both selfsimilar and non selfsimilar processes; and on the other hand a process can 
be selfsimilar without showing independence of the increments. A celebrated example 
of this second case is the so called fractional Brownian motion: this is a centered, 
Gaussian, ff-selfsimilar (for H £ [0, 1]) process B(t) with stationary increments, and 
covariance function 

U\2H , I \2H _\f _ s |2_ff 

E [B(t)B(s)} = 11 +l 1 2 1 l — , t, s > 0. 

It is apparent that this is nothing else than a generalization of the well known covari- 
ance function of the usual Brownian motion E [B(t)B(s)] = min(i, s) that is recovered 
when H = 1/2. Since B(t) is centered and Gaussian, this covariance function is all that 
is needed to define the process also if it is not Markovian. In fact a fractional Brow- 
nian motion coincides with the usual Brownian motion (and hence is Markovian with 
independent increments) only for H = 1/2, while for H ^ 1/2 it is non-Mar kovian, has 
correlated increments and for H > 1/2 shows long-range dependence. In other words 
an iJ-selfsimilar fractional Brownian motion with H ^ 1/2 is neither additive, nor 
Markovian: in fact it is not even a semimartingale, and hence few results of stochastic 
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calculus can be used. Prom another standpoint (see [7] p. 230) we can say that the 
selfsimilarity can have different origins: it can stem either from the length of the dis- 
tribution tails of independent increments, or from the correlation between short-tailed, 
Gaussian increments, and the two effects can also be mixed. For more information 
about the fractional Brownian motion see [Tl] and [12] 

3.3 Selfdecomposable laws and Ornstein— Uhlenbeck processes 

Selfdecomposable laws appear also in another important context: they are the most 
general class of stationary distributions of processes of the Ornstein- Uhlenbeck type. 
Take a Levy process Z(t) with generating triplet (a 2 , c) and logarithmic characteristic 
x(u), and for b > consider the stochastic differential equation (for simplicity we take 
T= 1) 

dX{t) = -bX{t) dt + dZ(t), X(0) = X P-a.s. (16) 

whose exact meaning is rather in its integral form 

X(t)=X -b [ X(s)ds + Z(t). 
Jo 

When Z(t) is a Wiener process the equation (fT6|) coincides with the stochastic differ- 
ential equation of an ordinary, Gaussian Ornstein-Uhlenbeck process; but the equa- 
tion (|16p keeps the meaning of a well behaved stochastic differential equation even if 
Z(t) is a generic, non Gaussian Levy process, and its solution 

X{t) = X e~ bt + f e b{ - s -^dZ(s). 
Jo 

will be called a process of the Ornstein-Uhlenbeck type. Of course to give a rigorous 
sense to this solution we should define our stochastic integrals for a generic Levy process 
Z(t): since all Levy processes are semimartingales (see [7J p. 255), this can certainly be 
done, but we will skip this point referring the reader to the existing literature (see [5] 
and |7j, or [13] for an extensive treatment). We will rather shift our attention to the 
possible existence of stationary distributions for a process of the Ornstein-Uhlenbeck 
type. In fact it is possible to show (see [5] p. 108, and [7J p. 485) that if 

/ log |x| /j,(dx) < +oo 

J\x\>l 

then the Ornstein-Uhlenbeck process X(t) solution of (|16p has a stationary distribution 
(f(u) = e^'") which is selfdecomposable with logarithmic characteristic 

r+oo 

V(u) = / x(ue- bt ) dt (17) 
Jo 

and generating triplet (/3 2 ,^, 7) where (3 2 = a 2 /2b, 7 = c/b, and - according to the 
equation ([9]) - the absolutely continuous Levy measure v has a density 

w(x) = ^ = — x / ^[ x > +0 °)}' ifx >°; 

|x| b\x\ \ fJ,{(— 00, x]}, ifx<0. 
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Conversely for every selfdecomposable law <p(u) there is a Levy process Z(t) such that 
ip(u) is the stationary law of the Ornstein-Uhlenbeck process driven by Z(t). Remark 
also that by the simple change of variable s = ue~ bt the relation (|17p takes the form 

^(«) = \l U — ds ' X(u) = bui//(u) (18) 

Jo s 

which is well suited to the inverse problem of finding the Levy noise of an Ornstein- 
Uhlenbeck process for a prescribed selfdecomposable stationary distribution. 



4 Examples 

4.1 Families of laws 

We will consider now several families of distributions (for further details see for exam- 
ple [10] and references quoted therein), all absolutely continuous, centered and sym- 
metric, with a space scale parameter a > which will of course span the types since 
the centering parameters always vanish: 

• the type of the Normal laws N{a) with density function and characteristic func- 
tion 



-x 2 /2a 2 



f(x) = ^-=- , <p{u) = e- a « / 2 

ay 2tt 



and with variance a 2 ; 

• the types (one for every A > 0) of the Variance-Gamma laws VQ(X, a) with 

(\x\/a) x -*K x _i(\x\/a) / i ^ 

a2 A " 1 r(A) V2vr \ 1 + a z u 2 J 

where K u (z) are the modified Bessel functions and T(z) is the Euler Gamma 
function [14] ; their variance 2Aa 2 is always finite; 

• the types (one for every A > 0) of the Student laws 7~(A, a) with density function 
and characteristic function 

, 1 ( a 2 \^ , , 2(a\u\)^K x/2 (a\u\) 



oB(i |) V« 2 + ^ 2 / ' 2V2r(A/2) 

where B(x,y) is the Euler Beta function p3j. Their variance is finite only for 
A > 2 and its value is a 2 /(A — 2). 

Important particular types within the Variance-Gamma and the Student families are 
respectively the Laplace (double exponential) laws C(a) = VQ{l,a) with density func- 
tion and characteristic function 

e -\x\/a j 
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and finite variance 2a 2 , and the Cauchy laws C(a) = 7"(l,a) with 

/(*) = — ^u)=e-^, 
air o z + x z 

and divergent variance. Finally we will also consider in the following another type of 
Student laws S(a) = T(3, a) with density function and characteristic function 



a7r \ a z + ar 



and finite variance a 2 . 

All the laws of our families are selfdecomposable (and hence infinitely divisible), but 
only N{a) and C(a) are types of a-stable distributions: more precisely jV(a) laws are 
2-stable, and C(a) are 1-st able. On the other hand the family VG(X, a) is closed under 
convolution, while T(A, a) is not. Of course this does not mean that the Variance- 
Gamma laws are stable since VQ(X,a) is not a unique type, and a convolution will 
mix different types with different A values. The infinitely divisible laws M(a), C(a), 
C(a) and <S(a) are of course entitled to their characteristic triplets (/3 2 , u, 7). Since they 
are all centered and symmetric we have 7 = for all of them. As for [3 2 it can be 
seen that for C(a), C(a) and S(a) we have j3 2 = (in fact, in terms of the Levy-Ito 
decomposition, they generate so-called pure jump processes), while for ]\[{a) we have 
f3 2 = a 2 . As for the Levy measures, on the other hand, we first of all have v = for 
the laws jV(a): from the point of view of the sample path properties this simply means 
that - at variance with the other three cases under present investigation - the Levy 
processes generated by Gaussian distributions never make jumps. The Levy measures 
of the other three cases are instead all absolutely continuous and have the following 
densities W(x) 



-\x\/a 



\X\ 

a 



irx 2 



7TX 2 



\X\ / \X\ \X\ \X\ \X\ 

1 I sin — ci cos — si — 

a \ a a a a 



for C(a) 
for C(a) 
for S(a) 



where the sine and the cosine integral functions for x > are |14j 

f +co sint , f + °° cost , 
si x = — / at , ci x = — / at . 

J X " J X " 

Examples of these three densities are plotted in Figure [TJ 



4.2 Convergence of consecutive sums 

To give examples of consecutive sums converging to our laws let us first of all recall 
what happens in the case of the Poisson laws V{X). Let B(n,p) represent the binomial 
laws for n independent trials of verification of an event occurring with probability p, 
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0.5 1 1.5 2 



Figure 1: Densities W(x) of the Levy measures for the selfdecomposable laws C(o) (black), 
C(o) (red) and S(a) (blue). To make the plots comparable we have chosen a = 1 for all the 
three densities, and to make the differences more visible we have plotted only the positive x-axis 
since the curves are exactly symmetric on the negative axis. 

and take the triangular array of Bernoulli 0-1 random variables ~ £>(1, A/n) with 
k = 1, . . . , n and n = 1,2,... Apparently they are uniformly, asymptotically negligible 
because for < e < 1 we have 

maxPf|xf ) |>e) = -^0. 
k I K ) n 

It is very well known that the consecutive sums are Binomial random variables, namely 

Sn = x[ n) + ... + X<p~B(n,£\ 

and that, according to the classical Poisson theorem, the limit law of these S n is V{X). 
Remark that, since in passing from n to n' the random variables change type, 

it will not be possible to put S n in the form of a normed sum as ([!]). The Poisson 
laws, however, are also limit laws in a still legitimate, but rather trivial sense due to 
the composition properties of the family 'P(A): take for instance a triangular array of 
Poisson random variables X^ 1 ' ~ V(X/n) with k = 1, . . . , n and n = 1,2,... which are 
again uniformly, asymptotically negligible because for < e < 1 we have 

max P {|4 n) | > e} = 1 - e~ x / n . 

Now - at variance with the previous example of the Bernoulli triangular array - at 
every step n we exactly have S n ~ ^(A) and hence, albeit in a trivial sense, the limit 
law again is V(X). Also in this case X% an d ^ belong to different types, so that 
it will be impossible to put S n in the form of a normed sum as ([!]). Of course both 
these examples show in what sense the Poisson laws are infinitely divisible but not 
selfdecomposable: they are limit laws of consecutive sums ([2]) of uniformly, asymptoti- 
cally negligible triangular arrays, but not of normed sums of independent random 
variables. 
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At the other end of the gamut of the infinitely divisible laws we find the 2-stable 
normal laws. To see in what sense they are limit laws take an arbitrary sequence X/. of 
centered, independent and identically distributed random variables with finite variance 
<t 2 and define the triangular array 

x (n) _ X k 

which always turns out to be uniformly, asymptotically negligible because of the Cheby- 
shev inequality: 

max p\\X ( k n) \ > e) = max P {\X k \ > ea^} < . 

The consecutive sums are now also normed sums of independent and identically dis- 
tributed random variables since 

1 n 

and according to the classical, normal Central Limit Theorem their limit law is AA(1). 
Also in this case, however, it is possible to exploit the composition properties of the 
normal type to find another, more trivial form of the consecutive sums: take the se- 
quence of normal independent and identically distributed random variables X k ~ AA(1) 
and define the triangular array xj^ = X k /yjn ~ N{l/n) which again apparently is 
uniformly, asymptotically negligible. Now the consecutive sums are also normed sums 
which - at variance with the previous example - are all normally distributed 

1 n 

and hence, in a trivial sense, the limit law is Finally let us remark that, since 

the normal laws are infinitely divisible, nothing will forbid them to be also limit laws 
of consecutive sums of triangular arrays that do not reduce to normed sums of inde- 
pendent and identically distributed random variables. However we will not elaborate 
here examples in this sense. 

We have introduced the trivial forms of the consecutive sums in the case of the 
Poisson and normal laws (see also the remarks at the end of the Section I2.1|) only 
because in our other subsequent examples this will be the unique explicit form available 
for our S n . Remark also that these trivial forms essentially derive from the fact that our 
laws are all infinitely divisible. In fact if (p is infinitely divisible, then also cp n = p l ^ n 
is a characteristic function, and of course <p™ = ip for every n. In general - with the 
exception of the stable laws - the ip n are not of the same type for different n, and 
hence the sums can not take the form of normed sums. That notwithstanding, in a 
trivial sense, every infinitely divisible law (p is the limit law of the consecutive sums 
of independent and identically distributed random variables all distributed according 
to ip n . What is less trivial, however, is to give an explicit form to the cumulative 
distribution function or density function of the component laws ip n : as the subsequent 
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examples will show this can be easily done only when we deal with families of laws 
closed under composition and decomposition. 

Let us consider first the Cauchy laws introduced in the previous Section I4.lt they 
are 1 -stable and, by taking advantage of the fact that the family C(a) is closed under 
composition and decomposition, it will be easy to show how they are limit laws of 
suitable sums of random variables. Take for instance a sequence of independent and 
identically distributed Cauchy random variables ~ C(o) and define the triangular 
array = X^/n ~ C(a/n). Since now there is no variance to speak about, to show 
that this sequence is uniformly, asymptotically negligible we can not use the Chebyshev 
inequality. If however 

11 x 

Fix) = — I — arctan — 
2 7r a 

is the common cumulative distribution function of the Xk, it is easy to see that the 
sequence is uniformly, asymptotically negligible because 

max P\\X { ™ ] \ >e} = max P{\X k \ > en} = 2 [1 - F(en)] 0. 

k I J k 

Now, as for the Gaussian case, the consecutive sums are also normed sums of indepen- 
dent and identically distributed random variables and are all distributed according to 
the Cauchy law C(a) 

1 n 

S n = X[ n) +... + lW = -Vl r C(a) 

so that, in a trivial sense, the limit law is C(a). What forbids here the convergence to the 
normal law is the fact that the variance is not finite, so that the normal Central Limit 
Theorem does not apply. At variance with the Poisson and Gaussian previous examples, 
however, we do not know non trivial forms of a Cauchy limit theorem embodying 
the stability of the Cauchy law. In other words we have neither explicit examples, 
nor general theorems characterizing the form of the normed sums of independent and 
identically distributed random variables whose laws converge to C(a), without being 
coincident with C(a) at every step n of the limiting process. 

This last remark holds also in the case of the Laplace selfdecomposable, but not 
stable laws C(a). In fact, since the VQ(X, a) family is closed under composition and 
decomposition, we can always take a triangular array X { ™ ] ~ VG(l/n, a) for k = 
1, ...,n and n = 1,2,..., and remark first that they are uniformly, asymptotically 
negligible by virtue of the Chebyshev inequality (the X^ have finite variance 
0), and then that for every n 

S n = X[ n) +... + XW- VG(1, a) = £{a) 

so that the limit law trivially is C(a). It must also be said that in this example the 
random variables of the triangular array change type with n so that the corresponding 
consecutive sums S n can not be recast in the form of normed sums of independent 
random variables. Since however the Laplace laws are not only infinitely divisible, but 
also selfdecomposable we would expect to find normed sums of independent (albeit not 
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identically distributed, because the Laplace laws are not stable) random variables whose 
laws converge to C(o). Unfortunately we do not have general theorems characterizing 
the needed sequences of independent random variables, and we can just show an ex- 
ample slightly more general than the previous one. Take for instance the triangular 
array with laws VQ ( p^rb^ri , a ) an d variances p^-T . They are uniformly, 



asymptotically negligible because from the Chebyshev inequality we have 

9 9 



{ 



max P < I X,. I > e > < max 



k e 2 k(2 + logn) e 2 (2 + logn) 
while, for a known property of the harmonic numbers, the consecutive sums are 



s n = x[ n) + ... + x^~vg f— i — 

1 » 2 + logn ^ k 



C(a) 



Now the sums S n are not trivially distributed according to C(a) at every n, but again 
they can not be put in the form of normed sums as they should since C(a) is selfde- 
composable. 

Finally similar remarks can be done for the selfdecomposable, but not stable Student 
laws S(a) introduced in the Section H~T1 but in this last case it is not even possible to give 
a simple form to the trivial consecutive sums because the Student family T(A, a) is not 
closed under composition and decomposition. In other words if ip is the characteristic 
function of a law S(a) we are sure that (p n = ip 1 ^ again is the characteristic function 
of a infinitely divisible law such that (p™ = <p for every n, but these component laws 
<p n no longer belong to the T(A, a) family as happens for the Variance-Gamma family, 
and in fact the form for their density function is rather complicated [10] . 

4.3 Stationary and selfsimilar processes 

Since all the laws of our examples are infinitely divisible we can use all of them to gen- 
erate the corresponding Levy processes by using (|12p to give the law of the increments 
on a time interval of width t. In particular we will explicitly do that for the laws M(a), 
C(a), C(a) and S(a). We then get as stationary increment characteristic functions <j>t{u) 
respectively 

e ~a tu /2T Wiener process from N(a) 

(l + a 2 u 2 ) t ^ T Laplace process from £(a) 

e ~at\u\/T Cauchy process from C(a) 

e~ a *H/ T (i + a \u\) t/T Student process from S(a) 

It is then apparent that the laws of the increments for the a-stable Wiener and Cauchy 
processes are simply M{a^J t/T) and C(at/T), while for the Laplace process the law of 
the increments is actually a Laplace law only for t = T, while in general it is a VQ(t/T, a) 
at other values of t. For our Student process, on the other hand, the situation is less 
simple because the Student family is not closed under convolution, and the increment 
law no longer is in T(X, a) for t ^ T. In this case it is not easy to find the actual 
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distribution from its Fourier transform (j>t(u), and only recently it has been suggested 
that the increments are distributed according to a mixture of other Student laws (for 
further details see [TO]). The Wiener and the Cauchy processes are fZ-selfsimilar with 
H = 1/2 and H = 1 respectively. This can also be seen by looking at the interplay 
between the two - spatial and temporal - scale parameters a and T. In fact in the 
Wiener and Cauchy processes these two scale parameters appear in two combinations 
- respectively a 2 /T and a/T - such that a change in the time units can always be 
compensated by a corresponding, suitable change in the space units; as a consequence 
the distribution of the process is left unchanged by these twin scale changes. This, on 
the other hand, would not be possible in the Laplace and Student processes since a and 
T no longer appear in such combinations. 

That notwithstanding we can achieve selfsimilarity in additive, non stationary pro- 
cesses produced by all our selfdecomposable laws. From (fl~5l) in fact we can give the 
characteristic function <f> s ^(u) of the increments in the interval [s,t] for our four types 
of law. First of all from the Normal type jV(a) we get 

namely 



+2H „2H 

X{t) - X(s) ~ M ( a 



2 ~ 8* 



These laws define processes which coincide with the usual Wiener process if and only 
if H = 1/2. For H ^ 1/2, on the other hand, our process is additive, iJ— self similar, 
and Gaussian with non stationary increments, and hence does not even coincide with 
a fractional Brownian motion which has stationary and correlated increments. In a 
similar way from the laws of the Cauchy type C(a) we get 



X(t)-X(s) ~C 



t H - s H 



and the process will coincide with the stationary (Levy) Cauchy process when H = 
1, while for H ^ 1 we have an additive, .ff-selfsimilar process with non stationary 
increments. From the Laplace type C{a) on the other hand we obtain an i?-selfsimilar 
(with H > 0), additive process when we take 



l + a 2 (£) u 2 Kt 



\ 2H 




) + 







1 



1 + a 2 (±) 2H 



so that, for s > 0, the law of the increment on an interval [s, t] actually is a mixture 
- with time-dependent probabilistic weights - of a law degenerate in x = and of a 
Laplace law: 



/s\2H 

X(i)-X( S )~(~) 5 + 



■ t 



c 



a 2 t 2H 

J<2H 



In other words this means that for s > there is always a non-zero probability that in 
[s, t] the process increment will simply vanish. Finally from the Student type S(a) we 
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get the additive, iJ— self similar, non stationary process with 

e -at*\u\/T« (l+s£| u |) 
4>s,t{u) = 7 v- , 

e-« H \ u \/ T " (l + f£|u|) 

so that X(t) ~ S (at H /T H \ while nothing simple enough can be said of the indepen- 
dent increment laws. 



4.4 Ornstein-Uhlenbeck stationary distributions 

All the Levy processes introduced in the Section [4.31 can now be used as driving noises 
of Ornstein-Uhlenbeck processes according to the discussion of the Section 13,31 Here 
we will only list the essential properties of the corresponding stationary distributions by 
analyzing their logarithmic characteristics (|18p . First of all, if b is the parameter of the 
process as in (I16p (remember that we took there T = 1 for simplicity), for the Wiener 
and the Cauchy driving noises we immediately have from (|18p that the logarithmic 
characteristics ip(u) of the stationary distributions are respectively 



a 2 u 2 
aid 



for a Wiener noise (usual Ornstein-Uhlenbeck process) 
for a Cauchy noise 



namely that the stationary distributions simply are N(a/ y2b) and C(a/b). In the case 
of the Wiener noise (namely in the case of the ordinary, Gaussian Ornstein-Uhlenbeck 
process) this means that the stationary distribution has a variance a 2 /2b while the 
Gaussian law generating the Wiener process had a variance a 2 . For the Cauchy noise, 
on the other hand, there is no variance to speak about. 

For the other two Ornstein-Uhlenbeck Levy noises (Laplace and Student) an im- 
portant role is played by the so called dilogarithm function |14^ \TE\ [TB] 




In fact a direct calculation of the integrals (1180 gives for the logarithmic characteristics 
tp(u) 

— Li2(— a 2 u 2 ) , for a Laplace noise 

_^Li _ Li2(— a\u\) , for a Student noise 

b b 

From the characteristic functions <p(u) = we can also calculate the stationary 

variances as —tp"(0) and we get a 2 jb and a 2 /2b respectively in the Laplace and in the 
Student case. Remark that - when they exist finite - the variances of the stationary 
distributions always are in the same relation with the variance of the law generating 
the noise: the stationary variance is the generating law variance divided by 26. The 
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Figure 2: Ornstcin-Uhlcnbeck process driven by a Laplace noise: characteristic functions 
(left) and density functions (right) of the stationary distribution (black lines), compared with 
the characteristic functions and density functions of the Laplace law C(a) generating the driving 
noise (red lines). Here a = \/\/2 and b — 1 so that both the variances (that of the stationary 
distribution, and that of C{a)) are equal to 1. 





Figure 3: Ornstein-Uhlenbeck process driven by a Student noise: characteristic functions 
(left) and density functions (right) of the stationary distribution (black lines), compared with 
the characteristic functions and density functions of the Student law S(a) generating the driving 
noise (red lines). Here a = b = 1 so that both the variances (that of the stationary distribution, 
and that of 5(a)) are equal to 1. 



form of the corresponding density functions is not known analytically, but it can be 
assessed by numerically calculating the inverse Fourier transforms of the characteristic 
functions: the results of these calculations for a couple of particular cases are shown in 
the Figures [2] and 

Finally, since the types of laws analyzed in this section are all selfdecomposable, by 
reversing the previous procedure we can also add a few remarks about the Ornstein- 
Uhlenbeck driving noises required to have Af(a), C(a), C(a) and S(a) as stationary 
distributions. In fact we have already said that for our two a-stable cases the stationary 
distributions are of the same type of the laws generating the driving noise, so that there 
is essentially nothing to add for the M(a) and C(a) stationary distributions. As for the 
other two cases on the other hand we will use the second equation (|18p to calculate the 
noise logarithmic characteristics x( u )'- 
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—2b tt-tt , for a Laplace C(a) stationary law 

1 + a z u z 

—b : — r , for a Student S(a) stationary law 

1 + a\u\ 

In both cases however the characteristic functions can not be elementarily inverted so 
that we do not have an explicit expression for the increment density functions of the 
the Levy noises that produce these two stationary Ornstein-Uhlenbeck distributions. 
We can only add a few remarks about the Laplace case: here the characteristic function 
<p(u) = e^ u ) does not vanish at the infinity since ip(±oo) = e~ 2b . As a consequence 
the Levy noise increment characteristic function can be better written as 

„2bt/(l+a 2 u 2 ) i 

\y(u)]* = e -2W« a « 3 /(l+a a « 3 ) = e -2W + (1 _ e -2bU £ £ 

so that the independent increments of an Ornstein-Uhlenbeck process with the Laplace 
type C(a) as stationary law are distributed according to a time-dependent mixture of 
two laws, one of which is degenerate in x = 0. As for the second law of this mixture, it 
has a density function given by 



I r+oa e 2bt/(l+a 2 u 2 ) _ j 

f(x,t) = - cos(nx) w du 

7T In e Zbt - 1 



but this integration can not be analytically performed. 



5 Conclusions 

Since many years selfsimilarity is a fashionable subject of investigation, in areas ranging 
from fractals to long-range interactions in complex systems: to have an idea just ask for 
the papers with the word "self similarity" either in their title or in their abstract present 
on arxiv. org and you will find 1 000 articles, and almost 200 of them only in the first 
six months of 2007. On the other hand this is a subject that has been approached 
from many different standpoints producing an unavoidable level of confusion [1] , while 
in fact it would be better discussed by placing the reader in the perspective of the 
general theory of the infinitely divisible (even non stable) processes. In the field of 
mathematical finance the use of non stable Levy processes is widespread, and several 
families of selfdecomposable laws and processes have been intensively studied in recent 
years: see for example the case of the Generalized Hyperbolic family [17 11181 [T9] . of the 
Student family pUGO], and of the Variance Gamma family (2TJ [22J [231 [23] . Considerable 
interest has also been elicited by the use of selfdecomposable laws in connection with the 
Ornstein-Uhlenbeck processes [3 [9], in particular for the stochastic volatility modelling. 
However, while in econophysics some non stable Levy laws are recognized as possible 
candidates for a consistent modelling of the underlying processes [25j[26], they remain 
less popular in the field of statistical mechanics and only recently their use has been 
proposed in connection with applications to the technology of accelerator beams [10^ 

In this paper we have tried to elucidate just a few points in the framework of the 
theory of stochastic processes. In particular we focused our attention on the relation 
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between on the one hand the selfsimilarity, and on the other the independence and the 
stationarity of the increments. This has led our inquiry toward the analysis of the laws 
of the process increments, and we have stressed the connection between the selfsimilar- 
ity of the process and the selfdecomposability of the increment laws. Selfdecomposable 
laws naturally arise in the study of the Central Limit Problem and of its solutions: 
in fact they are an intermediate (and more elusive) class of distributions between the 
more general infinitely divisible, and the more particular (and more popular) stable 
distributions. We found then that, in the case of selfdecomposable generating laws, 
both the stationarity of the increments and the selfsimilarity are always possible, but 
are not always present in the same process. On the other hand selfsimilarity can also 
be a property of (non Markovian) processes with non independent increments as in 
the case of the fractional Brownian motion. We finally stressed the connection between 
the selfdecomposability and the stationary laws of generalized Ornstein-Uhlenbeck pro- 
cesses with non Gaussian, Levy noises. All that has also be elucidated by means of a 
few particular examples, and some kind of application from physics to finance has also 
been pointed out. 
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